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O ■ Abstract 

Entanglement of spin and position variables produced by spatially inhomogeneous 
magnetic fields of Stern-Gerlach type acting on spinor Bose-Einstein condensates may 
. lead to interference effects at the level of one-boson densities. A model is worked 

out for these effects which is amenable to analytical calculation for gaussian shaped 
condensates. The resulting interference effects are sensitive to the spin polarization 
properties of the condensate. 

PACS number: 03.75.Fi 

The classic Stern-Gerlach experiment |I|] is currently understood as providing for the spin 
analysis of a beam of particles using the entanglement of spin and spatial degrees of freedom 
through the action of inhomogeneous magnetic fields and subsequent spatial analysis of 
the scattered flux|2[]. In this paper we explore effects resulting from this same type of 
entanglement when an essentially static spinor Bose-Einstein condensate with large coherence 
length is subjected to Stern-Gerlach field gradients. In the case of pure condensates that are 
sufficiently dilute so that mean field effects are a minor correction, the involved dynamics 
can be treated in a rather straightforward and simple way for appropriate, non-trivial field 
configurations using the ideas recently developed in ref. 

A simple but still rather general model state for a pure spinor (spin 1 for definiteness) 
condensate is one in which all bosons occupy one same single-boson spin-orbital u(f, a), 
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where a = — 1, 0, 1 is the spin variable. This state can be expressed in terms of the creation 
operator 



a J 

where ^> CT (f) is the spinor bosonic field operator, as 

This state is completely characterized by its one-boson density, given by 

p(f>;fV) = (¥ N) \tfjl(r')Mr)\^ {N) ) = Nu(r,a)u*(f',a'). 

The spin and position degrees of freedom are entangled whenever the spin-orbital u(r, a) 
does not have the product form u(r)x{ a )- A slightly more general form for a non-entangled 
one-body density is 

p(f,a;fV) ^u{r)u*{r')® p {s \a-a') = p {pos) (r;f') ® p (s) (a; a') (1) 

where the second factor, representing the spin density matrix, is allowed to be a mixed state. 
This information if of course not sufficient to craracterize the underlying many-boson state. 
Considering, as an example, a case in which p^{a;a') is diagonal, with eigenvalues N a and 
J2 a =-i N a = N, both the pure iV-boson state 

l(AW> = vw-i nw"' «i - / ^ «(o«ko (2) 

and the statistical mixture represented by the many-boson density 



E(4) \o)-£(o\M N (3) 



N 

lead to that same one boson density. 

The condensates recently produced using all-optical confinement techniques 0] may pos- 
sibly be associated with one-boson densities of type ([!]). This has been corroborated by 
applying a Stern-Gerlach magnetic field to the condensate and observing its splitting into 
three separate components. In general, the effect of Stern-Gerlach fields is to generate 
spin-position entanglement, by subjecting different spin components of the state to different 
boosts. One may note, moreover, that the spin components are defined with respect to a 
quantization axis associated to the field configuration itself, and are therefore open to ma- 
nipulation. In particular, a mid-course change in the quantization axis will imply that the 
position amplitude associated with each of the new spin components will be a superposition 
of differently boosted components of the initial state. This will give interference effects at 
the level of the one-body density. Alternatively, a final spin selective observation scheme 
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implemented in terms of a quantization axis different from that associated with the Stern- 
Gerlach field will also produce one-body interference effects. Both of these schemes will be 
developed in what follows. 

The mechanism responsible for the one-body interference can be illustrated schematically 
in terms of a completely polarized "travelling" one-boson amplitude 



pikr 



(r,a) = -j=®\a z = 0). 



If the field gradient is applied along the x-axis, rotation of the state changes the spin factor 
to (\a x = 1) — \o x = — 1))/V2. The Stern-Gerlach boosts then change the amplitude to (up 
to overall position independent phase factors) 



a) -> — = (e l ^ +iKX \a x = 1) - e t? ~ iKX \ 



where tin is the momentum associated with the boosts. The spin states \a x = ±1) can then 
be re-expressed in the z quantization clXIS cLS 

Wx = ±1) = \ (K = 1) ± V2\<r g = 0) + \a z = -1) 

so that the spatial density associated with the ^-component <j z = 1 (say) is, after the Stern- 
Gerlach entanglement 

p(f, a z = 1; r, a z = 1) = ^ sin 2 kx 

which shows interference fringes parallel to the z-axis. Note, in particular, that this is 
independent of the initial momentum k of the travelling amplitude. Eventually, an analogy 
may be discerned between this mechanism and that involved in Bragg diffraction experiments 
involving Bose-Einstein condensates 0, [|. More specifically, the coherent superposition of 
amplitudes with different momenta resulting from the action of the Stern-Gerlach field after 
subsequent redefinition of the relevant quantization axis parallels that achieved through the 
action of the Bragg pulses in the first ref.0. 

We now turn to the description of the dinamics of the spinor condensate under the field 
gradient. We assume that when the field is applied the condensate is dilute enough so that 
effects of the two-body interaction are negligible. The relevant dynamics is thus contained 
in a second-quantized Hamiltonian of the form 



-h 2 v 2 - - 

-<W + (a'\gS ■ B(f)\a) 
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The external magnetic field B(f) will be taken to be such that there is a suitable Cartesian 
frame in which it has the form, assumed to be valid throughout the spatial extention of the 
condensate (~ 100 (im), 
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B(f) — > B(x, z) = (£>o + h\£)z — b\xx. 

Typical values of the constants are Bq ~ 50 G and b\ ~ 2.5 G cmT x . 

As shown in ref [[J, the dynamics can be implemented in terms of a set of semi-classical 
propagators associated to eigenvalues of the projection of the spin operator onto the direction 
of the local magnetic field. In view of the orders of magnitude given above, this direction 
differs negligibly from the field z axis. Therefore, if this axis makes an angle 9 with the 
quantization axis used to express the spin density, in order to apply the semi-classsical 
propagators one has to consider the rotated spin density 

p^\m- m!) = £ < m (#)p (s V; *O«W(0) 

<7<J ! 

so that the complete initial density is given as 

p(r, m; r'm!) = p^ pos \f; r') £g> p^ s \m; m'). 

Propagation to time t of each m component is then performed by 

M \ i 



(x f ,Zf\U m (t)\x,z) w (— — -) exp ^5 c i as 



where the classical action 

M 2 2 1 (gmbi) 2 o 
5 c ias = - x ) + W - ■z) ] - 2^ m6lt ^ + 2 /) 24^"^ + ^ mjBot 

has been evaluated using m-dependent trajectories in the magnetic field. Note that m is a 
constant of motion along the trajectories due to the special choice of reference frame0. The 
evaluation of the above expression involved neglecting terms of the order (biAx/B ) 2 , with 
Ax of the order of the propagation range. Since the evolution in the y direction is free, it 
has been ommited from the above expressions. The propagated density at time t is therefore 
given by 



The spin inclusive observed spatial density after propagation in the Stern-Gerlach field is 
then given by the trace 

p(obs) ^ t) =Y^ p(r, m; r , m, t) . 

m 

In the case of an initial spin density which is diagonal in the magnetic field frame with weights 
3/5 : 1/5 : 1/5 for the components —1, 0, +1 and an initial gaussian spatial distribution one 
obtains the result shown in fig. |l|(a), for t = 15 ms and field gradient somewhat lower than 
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the typical value quoted above. The resulting distribution closely resembles that given in 
ref. 1. 

We next consider a mid-course change of the Stern-Gerlach field. After a first stage 
described by the above expressions and lasting for a time t\ a sudden rotation by an angle 
9b of B(x, z) is performed around the y axis. Strictly speaking, it is not important that the 
field values are maintained, so that the rotation can be achieved e.g. by superimposing an 
additional external field. The initial condition for the second stage of the evolution will be 
then, in the appropriate field frame 

p (2) (f,m;fW,ti) = dm im ^B)p(r,m 1 -,r , m 2 ,ti)d rn2ml (9 B ). 

miTO2 

The propagation to t > t% is handled analogously to the first stage, in terms of a propagator 
U m (t — ti), leading to a final density 

p (2) (f,m; f'm',t) = J d 3 n J d 3 r 2 {r\U m (t - ti)|ri)p (2) (ri, m; f 2 , rri, ti)(f 2 |E^*(t - h)\f '). 

The corresponding spin-inclusive final spatial density is now 

p (inc) ^ t) = P (2) ( f > ^ f ' m ' *) = ( 4 ) 

m 

= EE / d3r J ' ^{fjZ^, (t) If^joCr!, m l5 r 2 , rr^) (ra^^Ct) |r) 

m mim2 

where 

(r]U mm ,(t)\r') = J ^ 1 (r1f/ m (t-t 1 )|f 1 )^, m (^)(f 1 |f/ m/ (t 1 )|f') 

is the effective, spin non-diagonal propagator from time zero to time t. Note that, during 
the first stage, different m components are subjected to different boosts. After the rotation 
9b, the new m components will in general be associated with coherent superpositions of 
differently boosted spatial amplitudes, leading to interferece patterns in the observed spin 
inclusive one-boson density, as can be seen in fig. [I]. This interference will be absent in the 
special case in which the spin density is diagonal in the reference frame defined by the field 
acting in the first stage, due to the incoherent nature of the differently boosted components. 

Alternatively, we may let the system evolve as in the first stage for a time t, and then 
consider the resulting one-boson density spin-selected along a quantization axis differing 
from the field z axis by a rotation # se i around the y axis. In this case the spatial one-boson 
density selected for the component m of the spin in the last frame will be given by 

P {sel) (r,t) = d* mim (0sci)p(r,m 1 -,f,m 2 ,t)d m2m (9 se i). (5) 

which will also show interference effects for the same general reasons discussed in the previous 
case. In particular, in the special case in which the spin density is diagonal in the reference 
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frame defined by the field, no interference occurs in the one-boson density. An example of 
this case is shown in fig. |^. 

The initial (gaussian) spatial amplitudes leading to the patterns shown in Figs. 1 and 2 
have been taken as real, and therefore the initial velocity field vanishes. The regime in which 
the two-body interaction effects becomes negligible is in general reached after release of the 
mean-field energy, which generates a non- vanishing velocity field. The argument given above 
in order to illustrate the one-boson interference mechanism indicates that this initial velocity 
field does not change the interference pattern generated by the Stern-Gerlach boosts. This 
has in fact been verified in the calculations, by starting from narrower gaussians and letting 
them spread freely for the appropriate time before applying the magnetic field. 

In both cases, the typical fringe spacing is determined by the de Broglie wavelength 
associated with the relative velocities of the interfering amplitudes. This depends both on 
the value of the magnetic field gradient and on the duration of the first stage (in the first 
case above) or of the exposure time to the field (in the second case). For the F — 1 states 
of 85 Rb, fringes of the order of several microns correspond to times of the order of 1 ms or 
shorter for bi ~ 1 Gem -1 . 

It should be stressed that the above analysis has been carried out at the level of the 
one-boson density. To the extent that laboratory observation is sensitive to many-boson 
correlations 0], further interference effects are to be expected. Notably, in the case of the 
pure many-boson initial state (§) the Monte Carlo procedure of ref. @ will generate an in- 
terference pattern in the extended r, a space even when the initial Stern-Gerlach boosts are 
applied along the z-axis, which is however invisible at the one-boson level. No many-body 
interferences occur, under these conditions, in the case of the mixed many-boson state (|3|), 
due to the incoherece of the different spin components. If, however, the initial Stern-Gerlach 
boost is applied along a direction different from the x-axis, each one of the incoherent, 
completely polarized spin components of the mixed many-boson state will generate an in- 
terference pattern which seeps down to the one-body level. The restriction to the one-body 
level appears thus as a "filter" on the richer domain of many-body effects. 
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(a) t1=0. ms, t=15. ms 



(b) t1=.5 ms, t=1 5.5 ms 



Figure 1: Contour plots (sizes given in 
/im) of the spin inclusive one-boson den- 
sity (§) integrated over y for the shown 
values of ti and t. The initial one-boson 
spin density is diagonal with eigenvalue 
ratios 3 : 1 : 1 for spin ^-components 
—1,0,-1-1 respectively. The initial gaus- 
sian amplitude has b x = 25 fim and b VjZ = 
15 \im. The angles 9 and 9 b are 7r/2 and 
— 7r/2 respectively. Magnetic field param- 
eters: Bq = 50 G and b\ = 1.0 Gem -1 . 
Note that in (c) interference fringes are 
too narow for accurate representation. 
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Figure 2: Contour plots (sizes given in 

fim) of the m = +1 spin selected one- 2Q 

boson density @ integrated over y for the 3Q 

shown values of t. The initial one-boson 

spin density is a pure s z = — 1 state. The 

angles and 6 se \ are 7r/4 and n/2 respec- 60 

tively. Gaussian and field parameters are 50 

the same as in Fig. [l|. Note that in (b) 40 

interference fringes are too narow for ac- 30 

20 

curate representation. 
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